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We present a detailed study of the structural behavior and polarization reversal mechanism in 
phase III of KNO3, an unusual ferroelectric material in which the nitrate groups rotate during 
polarization reversal. This material was one of several studied in a previous work [O. Dieguez and 
D. Vanderbilt, Phys. Rev. Lett. 96, 056401 (2006)] where methods were described for computing 
curves of energy versus electric polarization. In the present work we extend and systematize the 
previous first- principles calculations on KNO3, and analyze in detail a two-parameter model in 
which the energy of the system is written as a low-order expansion in the polarization and the 
nitrate group orientation. We confirm that this model reproduces the first-principles results for 
KNO3 very well and construct its parameter-space phase diagram, describing regions where unusual 
triple-well potentials appear. We also present first-principles calculations of KNO3 under pressure, 
finding that its energy-versus-polarization curves change character by developing a first-derivative 
discontinuity at zero polarization. 

PACS numbers: PACS: 77.80.-e 81.05.Zx 
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I. INTRODUCTION 



Potassium nitrate (KNO3), also known as saltpetre, 
has long been used as an ingredient in explosives and pro- 
pellants, including 'black powder' and other early forms 
of gunpowder. At room temperature and pressure, KNO3 
crystallizes in an aragonite (Pnma) phase with four for- 
mula units. This phase is usually referred to as phase 
II, and its domain of existence is illustrated in the phase 
diagram of Fig. 1 in Ref . [J. When this material is heated 
at atmospheric pressure, a transition to phase I occurs. 
The structure of phase I seems to have been the subject 
of a controversy for years, with the current status being 
not much changed since the discussion given twenty years 
ago by Scott et al. in Ref. Q]. The proposed structure 
has global C3 symmetry but local Dq symmetry. Upon 
cooling at atmospheric pressure, phase I does not trans- 
form back to phase II; instead there is a narrow window 
of temperatures (aproximately from 113°C to 120°C) in 
which a ferroelectric phase, known as phase III, occurs. 
It is possible to broaden the range of temperatures in 
which phase III exists by applying hydrostatic pressure 
to the bulk system, or by growing it as a film. In this film 
form, KNO3 has been proposed as a promising material 
for random-access memory devices^ 

The unit cell of phase III KNO3 is a five-atom rhombo- 
hedral cell with RSm symmetry as depicted in Fig. Ufa). 
This structure can be viewed as resulting from the stack- 
ing, along the rhombohedral axis, of planes consisting 
alternately of K atoms and NO3 groups. However, it is 
found that a K plane is not equidistant from the NO3 
plane above it and the one below it. Instead, the system 
finds it favorable to adopt a polar structure in which the 
distance between a given K atom and its six O second- 
neighbors in the NO3 plane below it is reduced at the 
expense of increasing the distance between this same K 
atom and its three O first-neighbors in the NO3 plane 
above it [Fig. [IJd)]. Experimentally, the optimum stack- 



ing occurs when every K plane is roughly three times 
closer to the NO3 plane below it than to the one above.- 
Of course, there must also exist a polarization-reversed 
structure having the same energy but with the K plane 
closer to the NO3 plane above than to the one below. 
Moreover, the polarization reversal has to be accompa- 
nied by a NO3 rotation of 60° in order to preserve the 
symmetry. The resulting structure is shown in Figs, [ljc) 
and (f). A paraelectric configuration that is interme- 
diate between these two ferroelectric ones, having i?32 
symmetry, is shown in Fig. [TJb) and (e) ; in this case the 
K planes are equidistant from the NO3 planes above and 
below, and each K atom has six equidistant O first neigh- 
bors. This structure is not found to be stable in nature, 
but corresponds to the saddle point on the energy surface 
connecting the degenerate minima. If we take the nitrate 
orientation in the case of Fig. [TJa) to correspond to an 
angle of = 0°, then in the case of Fig. |ljb) we have 
= 30°, and in the case of Fig. [TJc) we have 6 = 60°. 

We presented preliminary results for ferroelectric 
KNO3 as an illustrative example in a previous study in 
which we developed a formalism for carrying out first- 
principles calculations at constant polarization^. We ob- 
tained curves of energy E versus polarization P for polar- 
izations parallel to the rhombohedral axis of KNO3 using 
not only the full first-principles method, but also an ear- 
lier approximate approach introduced by Sai, Rabe, and 
Vanderbilt (SRV)4 (The latter includes only the lattice- 
displacement responses to the electric field, while the 
method of Ref. y properly includes also the electronic 
contributions and is therefore exact.) We found that us- 
ing the SRV method led to an E{P) curve for KNO3 that 
was smooth, with the orientation of the NO3 group rotat- 
ing continuously from 6 = 0° to 9 — 60°. However, using 
the exact approach, we found a quite different E(P) curve 
having a cusp at zero polarization, reflecting the presence 
of a discontinuity in 6{P) at P = 0. A brief discussion 
of this was given in Ref. [H, where these features were 
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FIG. 1: (Color online) Top and side views of the fifteen-atom 
conventional hexagonal cell (outlined in black) for ferroelec- 
tric KNO3 in (a) the ground state with negative polarization, 
(b) a metastable paraelectric state, and (c) the ground state 
with positive polarization. Blue lines indicate the equiva- 
lent five-atom rhombohedral cells. Panels (d)-(f) show top 
views of the same structures as in (a)-(c), clarifying the near- 
neighbor environment of a K atom. Neighbors in the plane 
above and below are indicated with black and blue lines re- 
spectively; first K-O neighbors are drawn as dashed lines, 
while second K-0 neighbors are designated by dotted lines. 



explained through the introduction of a simple model in 
which the energy of the system was expanded to low or- 
der in the relevant degrees of fredom (the polarization P 
and the nitrate group rotation angle 0). 

The purpose of the present theoretical study is to in- 
vestigate the structural and dielectric properties of phase 
III of KNO3 in greater depth, and to provide a framework 
for understanding this and related systems. To do this, 
we combine extended first-principles calculations with a 
detailed analysis of the analytical model to be presented 
shortly. The picture that emerges is that of a system with 
a much richer variety of possible behaviors than can be 
found in conventional ferroelectrics. Some of these be- 
haviors are realized for KNO3 under external pressure, 
while others correspond to conditions that we have not 
reached in our first-principles calculations, but that could 
exist for related materials. 

The rest of this article is organized as follows. In 
Sect. Ill] we analyze the analytical model presented in 
Ref. y|. Section [Til] contains first-principles results for 
the ferroelectric phase of KNO3 under pressure. Finally, 
we summarize our results in Sect. IIV1 



II. ANALYTICAL MODEL 

In the model introduced in Ref. y, the energy may be 
written as 

E{P,6) = E + E* cos 69 

+A cos 129 + BP cos 39 + CP 2 (1) 

where E* serves as an energy scale for the definition of a 
dimensionless energy E = (E — Eq) / E* . We also define 
a polarization scale P* = y E* /C and a dimensionless 
polarization P = P/P*, and rewrite Eq. ([T]) in rescaled 
form as 

E{P, 9) = cos 69 + a cos 129 + 0P cos 36* + P 2 (2) 
where 

.4 



and 



(3- 



E* 



B 



VCE* 



(3) 



(4) 



are two dimensionless parameters determining the behav- 
ior of the model. 



A. Phase diagram 

We begin by analyzing in detail the possible behaviors 
that the model system of Eq. © can take as a and (3 are 
varied, thereby building up the phase diagram of Fig. [5J 
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FIG. 2: (Color online) Phase diagram of the model given by 
Eq. ([2]) including a paraelectric phase (i?32) and two ferro- 
electric phases (R2>m and R3). Each thick (thin) continuous 
line represents a first-order (second-order) phase transition. 
The broken lines enclose areas of the diagram mentioned in 
the text. The empty circle (square) corresponds to a fit of 
the exact (SRV) E{P) fixed-cell first-principles calculations 
for KN0 3 to Eqs. (HH); the filled symbols and the arrows 
indicate how the situation changes when stress in the cell is 
relaxed (zero pressure). 




FIG. 3: Energy versus polarization curves in area Ila (top 
curve), area lie (bottom curve), and on the boundary sepa- 
rating them (middle curve). The value of a is the same in all 
three cases; the system goes from paraelectric to ferroelectric 
with increasing (3. 



another equivalent discontinuity at a positive P = P' 
with 9 jumping from 60° to 60° — 9' , where P' > and 
a < 9' < 30°. The areas that arise in this way have 
been labeled by appending a letter (a, b, c, or d) to the 
Roman numeral used in the main classification. 



B. Triple- well potentials 



We will assume that (3 > 0, since the /? < case is equiva- 
lent to changing P into —P. We identify the local minima 
of E(P, 9) on the basis of its first and second derivatives, 
arriving at the following classification. First, there is a 
minimum at (P = 0, 6 = 30°) iff a < (8 - /3 2 )/32 (areas 
I and II in Fig. [J). Second, there are two equivalent min- 
ima at (P = -{3/2, 9 = 0°) and (P = +{3/2, 9 = 60°) iff 
a > ((3 2 - 8)/32 (areas II and III in Fig. \2§. And third, 
there are two equivalent minima with rotation angles sat- 
isfying cos6(9 = (/3 2 - 8)/(32a) and P = -(j3/2) cos30 
iff a > {{3 2 - 8)/32 and a > (8 - {3 2 )/32 (area IV in 
Fig. Note that local minima of the first and second 
kind coexist in area II, and we subdivide this according 
to whether the the global minimum is of the first kind 
(for {3 < 2V2) or of the second kind (for {3 > 2y/2). 
One is thus led to the phase diagram depicted in Fig. [2j 
showing three ground-state phases: a paraelectric i?32 
phase (area I and the left part of area II), a ferroelec- 
tric i?3 phase (area IV), and a ferroelectric Rim phase 
(the right part of area II and area III). The transition 
from phase i?32 to R3m is of first order, while the other 
transitions are of second order. 

We establish a further subclassification of the areas 
in the diagram by considering whether or not the cor- 
responding E(P) curves have continuous derivatives (or 
equivalently, whether or not 9{P) remains continuous). It 
can be shown that this depends only on a. For a > 1/4, 
there is a discontinuity in 9{P) at P = 0, with the angle 
jumping from 30° - 9 to 30° + 9 , where < 9 < 30°. 
For -1/20 < a < 1/4, 9(P) is continuous and E{P) is 
smooth. Finally, for a < —1/20 there is a discontinuity 
in 9{P) at P = —P' with 9 jumping from to 9' , and 



A very interesting area of the phase diagram of Figgis 
area II, in which a minimum at (P = 0, 9 = 30°) coexists 
with the pair of equivalent minima at 9 = 0° and 9 = 60°. 
It can be shown that if {3 < 2\/2 then the paraelectric 
minima has the lowest energy, while if (3 > 2\[2 the fer- 
roelectric minima are energetically favoured. Curves of 
E(P) for both cases, as well as for the {3 = 2^/2 case, 
are drawn in Fig. [3J They correspond to areas Ila and 
lie, in which the derivative for each of them is contin- 
uous, and illustrate clearly the first-order character of 
the transition from the i?32 phase to the R3m phase. 
(Analogous curves can be found at areas lib and lid, but 
in those cases there are two points at which the first- 
derivative of E{P) is not continuous.) These curves are 
similar to the ones that describe the free energy of a fer- 
roelectric material at temperatures in the vicinity of a 
first-order paraelectric- ferroelectric transition (see, e.g., 
Ref. d). However, we are not aware of any previous in- 
stance in which a first-principles mapping of energy ver- 
sus polarization has revealed this kind of triple-well be- 
havior at zero temperature. 



C. Fitting previous first-principles results 

The SRV results in Ref. d provide an example of a sys- 
tem showing a behavior that belongs in area Illb, while 
the exact results case falls into area Ilia. Fitting the first- 
principles results to Eq. |T]) and using the rescaling of 
Eq. ©, we find a = 0.105 and (3 = 4.06 for the SRV case, 
while instead a — 0.308 and (3 — 6.86 for the exact case. 
These points have been indicated in Fig. [2] with an empty 
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FIG. 4: Energy versus polarization, and NO3 rotation an- 
gle versus polarization, for the ferroelectric phase of KNO3 
without strain relaxation. The squares and circles represent 
SRV and exact first-principles calculations respectively (from 
Ref. 0), while the curves have been fitted to these data using 
Eqs. (QJ2J) as described in the text. 



1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 



- ° SRV 
: o EXACT 




-0.2 -0.1 0.0 0.1 0.2 

Polarization (C/m 2 ) 
-0.2 -0.1 0.0 0.1 0.2 



1 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 




■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ i ■ ■ ■ 



FIG. 5: Energy versus polarization, and NO3 rotation an- 
gle versus polarization, for the ferroelectric phase of KNO3 
with strain relaxation. The squares and circles represent SRV 
and exact first-principles calculations respectively, while the 
curves have been fitted to these data using Eqs. (UJ[2]) as de- 
scribed in the text. 



square and circle, respectively. The model curves that re- 
sult when these values are used, together with the points 
obtained from first-principles calculations, are plotted in 
Fig. |U We can see that in the region between the two 
minima the fit is excellent, indicating that our relatively 
simple model captures the important E(P) physics of the 
system. As the absolute magnitude of P grows too large, 
more powers of P are needed in the model to improve 
the fit. We have found that including terms up to fourth 
order works very well even for large values of P. 



III. FIRST-PRINCIPLES RESULTS 
A. Cell vectors relaxation 

The calculations in Ref. were carried out with the 
lattice vectors fixed to be those of the ground state, 
even when the polarization is varied away from its spon- 
taneous value. It is also possible to relax the cell 
vectors ^ we have now done that here, arriving at results 
valid for the more realistic situation of stress- free (i.e., 
zero-pressure) boundary conditions. As in our previous 
work;* we use density-functional theory 5 as implemented 
in the ABINIT code^ the local-density approximation 
(LDA) , 8 a plane- wave cutoff of 30 Ha, a reciprocal space 
grid with 6 inequivalent points, and Troullier-Martins 9 
pseudopotentials. 10 

The fitting parameters a and (3 resulting from these 
new calculations are shown in Fig. [2] as filled symbols, 



while the corresponding E(P) and 0(P) curves are shown 
in Fig. [SI These values are compared with the ones for the 
non-relaxed case in Table |TJ The SRV results are similar 
to the ones in the non-relaxed case, apart from a reduc- 
tion in the height of the double-well barrier that might 
have been expected in view of the additional degrees of 
freedom that are relaxed in the free-stress case. However, 
the exact results have changed qualitatively. The cusp at 
zero polarization and the corresponding discontinuity in 
the NO3 rotation angle have now disappeared, and the 
exact curves are quite similar to the SRV ones. We have 
also found that as the polarization crosses through zero, 
the cell volume expands appreciably. Specifically, the 
rhombohedral lattice parameter goes from 7.62 a.u. in the 
ferroelectric ground state to 8.00 a.u. in the paralectric 
case (a 5% increase) , with the rhombohedral angle going 
from 78.0° to 76.8°. This increase in lattice volume oc- 
curs so that the N-0 bond length can keep similar values 
in the ferroelectric i?3m ground state (2.353 a.u.) and in 



TABLE I: Values of the model parameters a and /3 depending 
on whether the cell vectors are relaxed or not, and whether 
we use the SRV method or the exact one. 



SRV Exact 

a a (3 

Fixed cell vectors 0.105 4.06 0.308 6.86 

Relaxed cell vectors 0.010 3.72 0.019 4.41 
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ally simpler than doing calculations at constant pressure, 
and we have confirmed that the anisotropy in the stress 
tensor obtained in this way remains quite small (~3 GPa 
at most), so that the rhombohedral angle would be ex- 
pected to change very little if the rhombohedral angle 
were relaxed. The results, shown in Fig. [6l indicate that 
when the lattice parameter is reduced by around 4.3%, 
a cusp reappears in the E(P) curve and the -R32 struc- 
ture becomes unstable to the R3 structure (9 ^ 30°.) at 
P = 0. The figure also shows that the N-0 bonds are 
not so tightly compressed as they would be at 9 = 30°. 
Therefore, we recover the Ilia and Illb behaviors de- 
scribed in Ref. d, but this time it is the pressure applied 
to the material that drives the transition between them. 



FIG. 6: (Color online) Zero-polarization first-principles calcu- 
lations of the value of the NO3 rotation angle and of the N-O 
distance in KNO3 as functions of the reduction in the lattice 
parameter. The empty symbols correspond to calculations 
in which the R32 structure is imposed, while the the filled 
symbols correspond to calculations in which this restriction is 
removed. All calculations were done using the exact method 
of Ref. 3 



the less closely-packed paraelectric i?32 state (2.351 a.u.). 
The reason for the appearance of the cusp at P=0 in the 
case of fixed lattice constants of 7.62 a.u. can be traced to 
the fact that the N-0 bond distance becomes uncomfort- 
ably short at 2.337 a.u. if the i?32 symmetry is enforced. 
Instead, the system prefers to lower its symmetry to i?3 
via a rotation of the nitrate groups (from 9 = 30° to 
36.1° or 23.9°, depending on whether P — > + or CP), 
with the N-0 distance attaining a slightly larger length 
of 2.338 a.u. 



B. Calculations under pressure 

The description provided in the previous paragraph 
seems to indicate that if the ferroelectric phase of KNO3 
were compressed, we should expect to see the E(P) cusp 
and 9{P) discontinuity at zero polarization reappearing 
at some value of applied pressure. To test this hypothesis, 
we have performed zero-polarization calculations for lat- 
tice parameters smaller than 8.00 a.u. while keeping the 
rhombohedral angle fixed at 76.8°. This is computation- 



IV. SUMMARY 

To summarize, we report a theoretical study of the 
properties of ferroelectric phase III of KNO3. We have 
carried out first-principles calculations that show that 
the polarization reversal mechanism in this material is 
accompanied by a rotation of 60° in the orientation of 
the nitrate groups. When the pressure exerted on the 
system is high enough, the corresponding double- well po- 
tential shows a distinctive cusp at zero polarization, and 
the mentioned rotation is discontinuous. These features 
are very well reproduced by a simple model in which the 
energy is expanded to low order in the polarization and 
the rotation angle of the nitrate groups. The phase di- 
agram that arises from this model contains regions that 
can be realized by varying the pressure on ferroelectric 
KNO3, as well as other regions that might turn out to be 
relevant for related materials in the future. One particu- 
larly interesting region of the phase diagram describes a 
system having a triple-well structure as a function of po- 
larization, a feature not yet encountered in the context of 
first-principles studies of ferroelectric energy landscapes. 
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